Quantum order-by-disorder near criticality and the secret of partial order in MnSi 
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The vicinity of quantum phase transitions has proven fertile ground in the search for new quantum 
phases. We propose a physically motivated and unifying description of phase reconstruction near 
metallic quantum-critical points using the idea of quantum order- by-disorder. Certain deformations 
of the Fermi surface associated with the onset of competing order enhance the phase space available 
for low-energy, particle-hole fluctuations and self-consistently lower the free energy. Applying the 
notion of quantum order-by-disorder to the itinerant helimagnet MnSi, we show that near to the 
quantum critical point, fluctuations lead to an increase of the spiral ordering wave vector and 
a reorientation away from the lattice favored directions. The magnetic ordering pattern in this 
fluctuation-driven phase is found to be in excellent agreement with the neutron scattering data in 
the partially ordered phase of MnSi. 



PACS numbers: 74.40. Kb, 74.40.-n,75.30.Kz, 75.70.Tj 

For a long time MnSi was considered a textbook ex- 
ample of a Fermi liquid developing ferromagnetic order 
at low temperatures where small spin-orbit (SO) cou- 
pling leads to a long-wavelength helical modulation of 
the magnetic order [TJ [2] . This view changed when ex- 
periments found a radical departure from conventional 
metallic physics upon suppression of magnetism under 
hydrostatic pressure. The occurrence of this behavior 
over a wide pressure range, not just close to the transi- 
tion [3 j, suggests a new phase of matter. A clue to its 
nature was provided by neutron scattering |4j, which re- 
vealed an unusual partially ordered phase, not apparent 
in susceptibility or resistivity. While the helimagnetic 
order is characterized by resolution-limited peaks corre- 
sponding to spiral wavevectors Q ~ 0.037 A -1 parallel to 
the [111] directions, in the partially ordered phase the 
scattering signal spreads diffusely over a sphere of ra- 
dius Q ~ 0.043A -1 , weakly favoring the [110] directions 
and anti- favoring the [111] and [100] directions. This is 
suggestive of spirals that fluctuate in direction, avoiding 
pinning by the lattice. 

In itinerant ferromagnets without SO coupling, it has 
been understood for some time that quantum fluctua- 
tions can drive the magnetic transition first order at low 
temperatures [5j [6] . The same physics has been revealed 
more recently in non-analytic corrections [7 to the Hertz- 
Millis theory [8] of ferromagnetic quantum criticality. 
This latter approach has been extended to systems with 
small SO coupling [9] , explaining the weak first order be- 
havior seen in early experiments on MnSi [10 j. Various 
itinerant ferromagnets display not only weak first-order 
behavior, but also the emergence of unusual phases in 
the vicinity of putative quantum critical points [TTJ [12] , 
suggesting that the avoidance of naked quantum critical 
points represents a generic principle [13 . In homoge- 
neous ferromagnets the quantum fluctuations that drive 
the transition first order also stabilize an inhomogeneous 



spiral phase, preempting the first-order transition [7| [14] . 

Formulation of these ideas in the form of quantum 
order-by-disorder [14] (or Coleman- Weinberg) provides 
an attractive physical interpretation. Moreover, it guides 
one to a simple route through the calculations. The cen- 
tral idea is that certain deformations of the Fermi surface 
enhance the phase space available for low-energy quan- 
tum fluctuations and so self-consistently lower the free 
energy. This mechanism is familiar in condensed mat- 
ter [T5] . Unusually, here it is driven by fermionic rather 
than bosonic fluctuations. The framework suggests di- 
rect connections to various experimental probes. Fermi- 
surface deformations and reconstructions associated with 
the onset of the competing order can be observed in 
photoemission and recently, the measurement of entropic 
landscapes has proven a revealing probe of phase recon- 
struction near quantum critical points [T2] . 

The utility of the quantum order-by-disorder approach 
is exemplified in its application to MnSi. We argue that 
at low temperatures near to the putative quantum criti- 
cal point, the spiral wavevector shifts away from its am- 
plitude and direction in the ordered phase in order to 
open up and benefit energetically from extra phase space 
for low-energy particle-hole excitations. This drives an 
instability towards a phase characterized by spirals that 
fluctuate about new directions, with an angular depen- 
dence as measured by neutron scattering [4 . 

The helimagnetic phase of MnSi is well described by 
a mean-field theory combining Stoner ferromagnetism 
with SO coupling resulting from the lack of inversion 
symmetry of the cubic B20 crystal structure. This 
can be captured phenomenologically by a Ginzburg- 
Landau free-energy density T = (VM) 2 + roM 2 + 
2DM • (V x M) + with D the Dzyaloshinskii- 
Moriya interaction. T is minimized by planar spirals 
M(r) = M [n x cos(Qr) + n y sin(Qr)] with propagation 
vector Q = Dn z (Fig. la). While T is invariant under 
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continuous rotations of the spiral, higher order spin-orbit 
terms lead to small corrections that reduce the symmetry 
to the point-group of the crystal. For MnSi, the leading 
anisotropy is given by pQ SJ 7 — A^(c^M^) 2 (i = a, 6, c 
labels crystal axes) which, for A < 0, leads to the exper- 
imentally observed pinning in the [111] direction [2]. 

The coefficients of the Landau function may be deter- 
mined microscopically by an expansion of the free energy 
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where U is the contact interaction between electrons, \i 
the chemical potential, are the mean-field electron 
dispersions in the presence of a planar spiral magnetiza- 
tion obtained from the microscopic Hamiltonian 



~\ h(k) -tr^c^Ck^, (2) 



with quadratic electron dispersion, contact repulsion and 
SO interaction h(k). c\. v and c\^ v denote fermion creation 
and annihilation operators, n rv the occupation operator 
in real space, and cr v y the vector of Pauli matrices. The 
SO coupling is a consequence of the lack of inversion sym- 
metry of the crystal structure of MnSi and, in the crys- 
tal axes basis, of the form h c = a[k a , fc^, k c ] + a[k a (k 2 — 
kl),k h (kl-k 2 a ),k c {k 2 a -kl)} [16]. Note that in the Hamil- 
tonian Q only the higher order SO coupling, a, breaks 
the rotational symmetry and is therefore responsible for 
the directional dependencies of both the helimagnet and 
the partially ordered phase. 

The mean-field decoupling of the interaction takes 
the form H v « U^ r ,uy M ( T ) ' ^vyc\ v c YV > = 
UM / k ( c k + Q/2t Ck -Q/ 2 i n - c 0- We have chosen the 
propagation vector Q = Qn z = Q (sin 6 cos (j)n a + 
sin^sin^n^ + cos#n c ) as quantization axis and so ro- 
tate the spin-orbit interaction to the spiral basis; h(k) = 
R^h c (R^k) = ak + af(k, <j>,0) where R0,#n c = n z 
(see Fig. lb). The components of f are lengthy third- 
order polynomials in k x: k y , k z with coefficients that de- 
pend on the direction of Q. To leading order (linear cou- 
pling to Q, corresponding to the Dzyaloshinskii-Moriya 
interaction) only h z contributes, yielding dispersions 
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where QT X = [f z (k + §) - f z (k - §)]/2 and Q 2 r 2 = 
[/*(k+§)+/*(k-§)]/2. The isotropic SO coupling leads 
to the Dzyaloshinskii-Moriya interaction with D = a. As 
shown later, the leading, mean-field anisotropy is indeed 
of the form ST = A EiO^M;) 2 with A - a 2 . 




n z = Q/Q (b) 




FIG. 1: (color online) (a) Planar spiral M(r) = 
M [n x cos(Qr) + n y sin(Qr)] with propagation vector Q = 
Qn z . The magnetic moments of size M are confined to a 
plane spanned by n^n^ _L Q. (b) Due to higher order SO 
terms in the cubic B20 crystal structure of MnSi, the free 
energy depends on the direction of Q = Q(sin# cos (j)\\ a + 
sin sin + cos#n c ) with respect to the crystal axes n a , 
and n c . 



In order to analyze the stabilization of competing 
ground states by quantum fluctuations, we include the 
leading fluctuation corrections which in self-consistent 
second order perturbation theory are given by [14] 



F fl = -2U 2 J2 n(e ^ )n(ek2)[n(e ^ )+n(e ^)] 

ki+k 2 =k 3 +k 4 6k i ^ 6k 2 _ 6 k 3 _ 6 k 4 

(4) 

with n(e) = (e^ e ^>I T + 1) 1 the Fermi function. 

Since the anisotropic SO coupling, a, is an order of 
magnitude smaller than a, it has a negligible effect upon 
the location of the phase boundaries and the pitch Q [17] . 
Therefore, we calculate the magnetic phase diagram for 
a = and then analyze the directional dependencies. For 
a = 0, the mean-field dispersion is identical to that with- 
out SO coupling if we replace Q by Q = Q — a. Therefore 
the phase diagram is the same as that of the homoge- 
neous systems with the appropriate shift of Q. The ho- 
mogeneous system shows a ferromagnet to paramagnetic 
transition with a fluctuation induced spiral phase form- 
ing around the putative quantum critical point. In the 
presence of SO coupling, a helimagnet to paramagnetic 
transition is obtained with a small region of fluctuation 
induced partial order around the critical point. 

Restricting our consideration to planar spiral configu- 
rations, we seek an expansion of the fluctuation-corrected 
free energy of the form F = (a + a 2 Q 2 + a 4 Q 4 ) M 2 + (&o + 
b2Q 2 )M 4 + coM 6 , where the coefficients ao, &2, e ^ c - are 
functions of J7, T and \i. This expansion is controlled 
around a finite temperature tricritical point (that will be 
revealed below) at which M = 0. 

Determination of the phase diagram is simplified by 
relationships between the expansion coefficients: At low 
temperatures, derivatives of the Fermi functions are 
sharply peaked at the Fermi surface leading to the pro- 
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FIG. 2: (color online) (a) Magnetic phase diagram for weak 
SO interaction a as a function of inverse interaction u -1 = 
fi/U and temperature t = T/ji. For strong interactions, the 
system develops helimagnetic order of pitch Qo = ot. The 
transition to the paramagnet is continuous. For smaller u, 
quantum fluctuations stabilize a different spiral ground state 
which emerges from the tricritical point shown in red. The 
higher-order SO term a, which is responsible for the direc- 
tional dependencies, has a negligible effect upon the location 
of the phase boundaries and the pitch Q [T7j . (b) Magnetiza- 
tion M and relative pitch Q = Q — Qo as a function of u~ x 
for t = 0.05, 0.1, 0.15, 0.2. At the transition between the heli- 
magnet and the partially ordered phase, Q ~ (u -1 — Uc 1 ) 1 ^ 2 
while M is continuous. The transition to the paramagnet is 
characterized by a first-order jump in M. 



portionalities a2 ~ 2(fc^)6o, a a ~ 3(kj)co and 62 ~ 
3(fc?)co (k z — k z /k) with angular averages (k 2n ) = 
l/(2n + 1). The remaining coefficients ao, bo and Co for 
Q = can be computed numerically over the whole tem- 
perature range and analytically at low temperatures [IS] . 



In the limit T — >> 0, we obtain a™* 

4!(2tt) 2 ^ anCL C 



2^2 2 Lmf 
(2tt) 2 M ' ^0 



cients and clq 



= 4 .gf^)2 ^ 6 f° r the mean- field coeffi- 
A(l + 2 In 2)^ 4 and b% ~ A(l + In t)^ 6 
for the fluctuation corrections, with dimensionless units 
u = U//jL and t = T/ja and where A = ^ 2 ^) 6 • 

The phase diagram is shown in Fig. 2 as a function 
of t and u — we expect this to reproduce the topology 
of the temperature-pressure phase diagram. A second 
order phase transition between the helimagnet and the 
paramagnet occurs along the curve ao = 0. A crucial 
factor determining the phase diagram is that quantum 
fluctuations make a In t contribution to bo [19] , which be- 
comes negative below t* = 0.24 as a result. On its own, 
this would imply that the transition turns first order. In 
combination with the concurrent change in sign of a 2 ~ 
26o/3, this signifies the onset of fluctuation-driven re- 
orientation of the spiral order and an increase in the spi- 
ral pitch. The transition between the fluctuation-driven 
spiral and the paramagnet along the curve aoco = jj^&o * s 
first order and characterized by a jump in magnetization 
(Fig. 2b). The transition between the two spiral phases 
(ao = line at t < t c ) is a Lifshitz line where Q increases 
continuously. Deviations from a k 2 dispersion render the 
transition weakly first order [18] . Experimentally, a small 
increase in Q is indeed observed [4]. 



The dependence of the free energy upon the direction 
of Q results from the anisotropic higher order SO cou- 
pling a and the corresponding deformations T 1 / 2 of the 
electron dispersion. At low T, the resulting free-energy 
contributions are given by 8F\ = ^-a 2 ao(Tl)Q 2 M 2 and 
SF 2 = 2a 2 bo(T 2 )Q 4 M 2 , where (...) denotes an angu- 
lar average over the direction of k. 5F\ is dominant 
in the ordered helimagnet and precisely of the form 
A J d 3 r ^2 a (d a M a ) 2 pQ when evaluated for planar spirals, 



where A 
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aod . It is minimized for spirals point- 



ing along the [111] directions. Although SF 2 is of higher 
order in the SO coupling a ~ Q it dominates the direc- 
tional dependence in the partially ordered phase since it 
is proportional to bo which diverges logarithmically at 
low T. The angular average yields (r?,) = 0) with 

0(0,0) = — sin(20) 2 [sin(0) -3sin(3<9)] 2 
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-cos(20) 2 sin(20) 2 



(5) 



which has maxima in the [110] directions and vanishes 
along [111] and [100]. Evaluating SF = SF 1 + SF 2 for 
different spiral orientations and setting 6F[m] = 5-Fjiio] 
we obtain the condition ao = |fro<^ 2 for the directional 
change. Up to order a to which we obtained the phase 
diagram, this coincides with the phase boundary between 
the spiral phases given by ao = 0. 

Going beyond the Ginzburg-Landau analysis and al- 
lowing for thermal fluctuations, the system forms do- 
mains in which the spirals are oriented around one of the 
principle directions. Within these domains, the system 
forms a statistical mechanical ensemble of spirals with 
fluctuating and spatially modulated orientation about 
this principle direction. Since the partially ordered phase 
is driven by fluctuations, the angular distribution is ex- 
pected to be larger and the timescale of fluctuations to 
be shorter than in the helimagnetic phase. 

Allowing for slow spatial modulation of the spiral ori- 
entation, we find a contribution to the free energy den- 
sity given by F = §M 2 |6 |(VQ) 2 . Evaluating the fluc- 
tuations in orientation using this and the previously cal- 
culated anisotropy, we obtain the same form as using 
the Boltzmann weight exp[— V G s 5F 2 /T] with V e s(T) ~ 
(Mc 2 a- 2 |6 |" 3/2 T- 1 / 2 )exp[a 2 M 2 |6 | 3 /(cgT)]. At very 
low temperatures, V G s —> 00 and fluctuations in orienta- 
tion become small. In Fig. 3, this weight is shown along a 
great circle through the high-symmetry points and as in- 
tensity plot over the momentum sphere documenting the 
excellent agreement with the neutron-scattering data [4] . 
Since 5F\ and SF 2 are proportional to M 2 , and V e fi is pro- 
portional to M, all fall towards zero along the helimagnet 
to paramagnet transition and are significantly reduced in 
the partially ordered phase (Fig. 3b). We anticipate neu- 
tron scattering peaks that are much more diffuse in angle 
in both these regions of the phase diagram. 

The timescales of fluctuations in the direction of Q 
are harder to calculate. Presumably they can be ap- 
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FIG. 3: (color online) (a) Boltzmann weight I(<fi,9) = 
Io exp[— V e ffSF2/T] = I exp[cg(0, 0)] of the term SF 2 which 
dominates the directional dependenve in the partially ordered 
phase. The overall pref actor c enters as the only free fitting 
parameter. Fluctuations favor spirals along [110] and equiv- 
alent directions while almost no intensity is found along the 
[111] directions of the helimagnet. (b) /((/>,#) with c = 4 
along a great circle connecting the high-symmetry directions 
in quantitative comparison with the experimental data pQ. 

proximated by a Fermi Golden Rule calculation of the 
rate for a two-step process where the electron fluid ab- 
sorbs a low-energy phonon whilst producing a particle- 
hole pair, followed by recombination of the particle-hole 
pair in a state with rotated Q and emission of a phonon 
with shifted energy. The enhanced density of low-energy 
particle-hole states that drive the formation of the par- 
tially ordered phase would then naturally lead to faster 
angular fluctuations. This expectation is consistent with 
the fact that the partially ordered phase is not visible in 
NMR - a probe that averages the spin configuration over 
a longer timescale than neutron scattering. 

Discussion and Conclusions: Quantum order-by- 
disorder provides a natural explanation of the partially 
ordered phase of MnSi. In particular, the angular distri- 
bution identified from neutron scattering is reproduced; 
the leading anisotropy of the SO coupling responsible 
for the [111] ordering in the helimagnetic phase auto- 
matically leads to [111] (and [100]) directions being anti- 
favored by fluctuations. This had proven problematic in 
previous analyses — it was the motivation for an alterna- 
tive suggestion of skyrmion-like ordering [20] . Skyrmion- 
like spin crystals (formed by the superposition of spirals 
with axes pointing in different directions) seem neces- 
sary to understand features in Hall effect seen at finite 
magnetic field [2Tj and are compatible with our analy- 
sis. However, they are unnecessary to describe the phase 
boundaries and anisotropy of the partially ordered phase. 

Several measurements would distinguish between spi- 
rals and skyrmions: i. Annealing in a magnetic field 
applied along one of the principle directions would align 
the spiral domains. Neutron scattering would show an 
intensity imbalance of the scattering peaks - an effect 
not expected for skyrmions; ii. Annealed samples should 



also show the anisotropic transport found upon phase re- 
construction near to other quantum critical points [T2] ; 
iii. The Lifshitz-like transition between helimagnet and 
partial order is expected to show a diffuse spreading of 
scattering over the whole sphere along the transition line. 

The unusual, non-Fermi liquid T 3 / 2 resistivity [3] re- 
mains puzzling. It is suspected that the effects of disor- 
der upon the partially ordered phase provide the key to 
understanding it [23 . The consequences of this within 
the order-by-disorder approach is an intriguing avenue 
for further study. 

The notion that fluctuations in the direction of the spi- 
ral order might drive the new physics in MnSi has been 
suggested by several authors [9j [TTl [22] working with 
extensions of Hertz-Millis theory. Indeed, the blue fog 
phase of Ref. [22 is similar in broad concept to our pic- 
ture for the partially ordered phase. The complementary 
perspective afforded by the quantum order-by-disorder 
approach — focussing upon particle-hole fluctuations in 
the presence of spiral order rather than fluctuations of 
the order parameter — allows a more tractable calcu- 
lations and leads to some qualitatively different results. 
We anticipate a similar utility in the broader context of 
phase-formation near to quantum critical points. 
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